A classical bounce: constraints and consequences 
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We perform a detailed investigation of the simplest possible cosmological model in whicii a bounce 
can occur, namely that where the dynamics is led by a simple massive scalar field in a general self- 
interacting potential and a background spacetime with positively curved spatial sections. By means 
of a phase space analysis, we give the conditions under which an initially contracting phase can 
be followed by a bounce and an inflationary phase lasting long enough (i.e., at least 60-70 e-folds) 
to suppress spatial curvature in today's observable universe. We find that, quite generically, this 
realization requires some amount of fine-tuning of the initial conditions. We study the effect of this 
background evolution on scalar perturbations by propagating an initial power-law power spectrum 
through the contracting phase, the bounce and the inflationary phase. We find that it is drastically 
modified, both spectrally (fc— mode mixing) and in amplitude. It also acquires, at leading order, 
an oscillatory component, which, once evolved through the radiation and matter dominated eras, 
happens to be compatible with observational data. 

PACS numbers: 



I. INTRODUCTION 

The inflationary paradigm is now part of the standard 
cosmological model (see Ref. [1 for a recent review with 
updated references). So much so that the very idea of 
challenging its main hypothesis might appear hopeless 
in view of its impressive success in explaining the oth- 
erwise mysterious observations that the Universe is flat 
and seemingly free of remnants such as monopoles and 
in providing a mechanism that exponentially damps any 
initial anisotropy [2], as well as considerably alleviates 
the homogeneity problem, i.e., instead of demanding ho- 
mogeneity over an almost infinite spacelike hypersurface, 
it reduces to finding a Planck size region over which the 
infiaton field is homogeneous. In the usual (chaotic) sce- 
nario [3J Sj, this is almost certain to occur, taking into 
account the infinite amount of possible initial conditions. 
It then suffices to argue that we happen to live in what 
has become of this initial region (i.e., applying some sort 
of anthropic principle). Moreover, specific examples have 
shown that an initially inhomogeneous universe, of the 
Tolman-Bondi type for instance [5|, also shows a trend 
to homogeneize under the action of inflation. Therefore, 
and even though these arguments may not be entirely 
convincing in the most general case, it is plausible to ar- 
gue, somehow, that inflation provides a dynamical means 
to drive any "weird" universe into one satisfying the cos- 
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mological principle. Finally, inflation also predicts the 
spectrum of primordial density fluctuations thanks to 
which large-scale structures formed; this predicted spec- 
trum can easily be made to fit all the known data. 

Why, then, would one insist in finding an alternative? 

First of all, having a serious competitor usually boosts 
the understanding and the predictability of the defender. 
Although it is certainly of interest to provide new tests 
for the inflationary paradigm to confront, it is of crucial 
importance to provide a priori alternatives were the infla- 
tionary predictions unable to match future observations. 
In order to do so, one naturally turns to high energy ex- 
tensions of the standard theory. Nowadays, this mostly 
means looking at the cosmological consequences of mod- 
els that can be implemented in string theory (see, e.g., 
[6J) in a satisfactory way. 

Moreover, inflation itself is not free from difficulties 
(see for instance the discussion by R. Brandenberger in 
Ref. [1]). Aside from assuming an otherwise never ob- 
served scalar field to lead the dynamics of the universe, 
scalar field whose potential is subject to some amount of 
fine-tuning (arguably at the 10^^^ level), it usually im- 
plies the existence of a singularity in the far past, may 
face a trans-Planckian problem, and demands that quasi- 
classical General Relativity (GR) be valid up to energies 
of the order 10~'^mpj, with m^^^ = G„ the Planck mass 
(G„ being Newton constant). Of these problems, none is 
by itself sufficient to reject the paradigm, but the combi- 
nation might provide a good motivation either to use a 
different approach or to complement inflation somehow. 

The pre Big-Bang scenario (see the review [7]) was 
the first to set the issue along these lines, followed, after 
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quite some time, by an application of the brane idea in 
e.g., Ref. [8] and in the so-called "ekpyrotic" model (see 
for instance [51 HH])- All these models^, written in the 
Einstein frame, contain a contracting phase preceding 
the ongoing expanding one. The reversal between these 
phases yields a bounce in the scale factor [TS j [TBI [T7| [T8] . 
with an impact on cosmological observations for which 
general properties have not yet been obtained. Note that 
roughly at the same time as these models were developed, 
discussions based purely on GR had concluded that no 
bounce could be achieved with a perfect fluid source [T^ . 

Among previous works on bouncing cosmologies, those 
that stay within the framework of 4-dimensional GR pOl 
121] have focused solely on the bounce-producing mech- 
anism. With a single dynamical degree of freedom, a 
scalar field say, this requires closed spatial sections in or- 
der to satisfy the null energy condition all along. Other 
choices are possibles. Achieving a bounce in 4 dimen- 
sions with flat or even open spatial sections requires one 
to violate the null energy condition and the presence of 
at least two kinds of fluids, one of which must have neg- 
ative energy density [21 [Ml 123 As another 
alternative, a scalar field with a nonconventional kinetic 
term can be used, the so-called "K— bounce" [27]. There 
also exists phenomenological descriptions which aim to 
classify the possible perturbation evolution through the 
bounce [^1^. 

Having managed a background bounce, one needs to 
describe the evolution of scalar and tensor perturba- 
tions in this background for comparison with observa- 
tions (note also the existence of other possible tests of 
bouncing scenarios, e.g., Ref. [30]). Propagating pertur- 
bations through the bounce itself is very intricate [31] [32] 
as one wants perturbations to remain small at all times. 
Unfortunately, Einstein gravity is often spoiled by the 
necessarily weird conditions to be imposed in order for 
a bounce to actually take place. In particular, having 
a negative energy fluid clearly leads to instabilities in 
the long run; this is especially true for models having 
an infinitely long contracting phase. This can however 
be handled by assuming the bounce-making negative en- 
ergy component to be merely phenomenological, acting 
only for a short time, the bounce duration, in general 
assumed comparable with the Planck scale. In any case, 
there are no generic properties that can be derived in 
a model-independent way as in slow-roll infiation [53] . 
Very often, though, it is found that a reasonable model 
connects two almost de Sitter phases. 

The presence of an expanding de Sitter phase after 
the bounce naturally leads to the question of whether 
one can implement an inflationary epoch connected to 
a bounce phase, in a mixed situation somehow reminis- 



cent of the topological defects [53] versus inflation chal- 
lenge: as topological defects could not be identified as 
the sole source of primordial perturbations, in particular 
because of their incoherent spectrum, it was suggested 
to switch the question from "infiation or defects" to "in- 
fiation and defects or infiation only" [35], especially in 
view of the fact that most reasonable grand unified theo- 
ries [3n] are expected to produce such defects at the end 
of infiation [57]. The idea of the present work relies on 
the same point of view: if a purely bouncing alternative 
is found to be unrealistic, why not consider a situation 
in which infiation follows a bounce? 

The purpose of this work is thus to study a class of 
models implementing this idea. We assume a bouncing 
phase, obtained by a simple scalar field evolving in a 
universe with closed spatial section, followed by an infia- 
tionary epoch. This second era then provides the stan- 
dard solution for some usual cosmological issues (fiatness 
and homogeneity as well as reheating), while the bounce 
permits the model to avoid a primordial singularity and 
provides an infinite horizon. The hope is then to find 
out whether any effect could somehow be present in the 
primordial spectrum, leaving some imprint to be tested 
against observations. 

In the following section, § |llj we discuss the back- 
ground model and recall the necessary conditions for the 
occurence of a bounce; this section is heavily based on 
Ref. [20,. We also implement constraints that were not 
considered in this reference and characterize the poten- 
tial through which perturbations propagate (this is de- 
tailed in § IV I. § III is then devoted to describing how a 



bounce followed by infiation can take place in our frame- 
work by means of a phase space analysis. This analysis 
indicates that, as expected, some amount of fine-tuning 
is demanded. In § IV we obtain the potential through 



^ We do not here consider the other category of models based on 
string gas cosmology 1121 1131 114| . as they do not contain a 
bounce phase. 



which perturbations propagate and describe the propa- 
gation of scalar perturbations in this background. We 
derive the resulting spectrum, which turns out to be the 
product of an oscillatory component and an almost fiat 
power-law component at leading order, and find that the 
curvature perturbation Cbst (as first defined by Bardeen, 
Steinhardt and Turner [38 J is not conserved through the 
bounce. On the one hand, we can conclude that, based 
on the current belief that the spectrum of primordial per- 
turbations at horizon exit should be scale-invariant and 
almost fiat with possibly superimposed features at higher 
order, this class of models can most likely be generically 
ruled out unless a justification is found under which the 
wavelength of the oscillations is sufficiently long so as to 
make them unnoticeable. On the other hand, because the 
cosmic microwave background (CMB) multipoles Ce are 
roughly given by the convolution of a spherical Bessel 
function with the initial power spectrum we find that 
evolving this spectrum through the radiation and mat- 
ter dominated eras damps the oscillations and results in 
a spectrum of CMB multipoles not incompatible with 
the Wilkinson Microwave Anisotropy Probe (WMAP) 
data |39j|4Qj. 
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II. BACKGROUND BOUNCING MODEL 

Our starting point consists of GR with a scalar field, 
i.e., we assume the dynamics to derive from the action 



S 



R 



V{ip) 



(1) 



with — |7rG„ the Planck length, (G„ being New- 
ton constant) and V {(p) the self-interaction potential 
for the scalar field (/?, which is left unspecified for the 
time being. In what follows, for notational convenience, 
we shall assume the natural system of units in which 
SttG^ — 1. Assuming homogeneity and isotropy, the 
background which solves the equations of motion derived 
from Eq. ([ij takes the Friedmann-Lemaitre-Robertson- 
Walker form and reads 



dt^ 



dr2 



1 -/Cr2 



'de^ 



r^sin^ed(t)' 



(2) 

where ait) is the scale factor, and where the spatial cur- 
vature /C > can be normalized to unity. The scalar 
field (p can be understood as either fundamental or phe- 
nomenological, but in both cases, its energy density p^p 
and pressure are given by 

• 2 '2 

P^^\ + V{v), p^^^^Viip). (3) 

In Eq. (|3| and in the following, a dot denotes a deriva- 
tive with respect to cosmic time t. In the forthcoming 
calculations, we also use the conformal time rj defined 
through 



dt = a{rj)dri, 



(4) 



and unless specified otherwise, derivatives w.r.t. rj will 
be denoted by a prime, so that for an arbitrary function 
of time /, one has /' = af. 



Einstein's field equations relate the time evolution of 
the scale factor a{t) to the stress tensor for ip (i.e., to the 
pressure and the energy density p^) through 



H' = 



1 /I 



3 V2 
1 



2 



(5) 
(6) 



where H = a/a is the Hubble expansion rate, while the 
(redundant) dynamical Klein-Gordon equation for p is 
given by 

ip + 3Hip + V^^0. (7) 
Combining Eqs. Q and Q yields the following relations 



Pv^+P^ = 2 -ijj , (8) 

P^ + 3p^ = -6 (ij + , (9) 
with which the energy conditions can be rephrased in 
terms of the time behaviour of the Hubble rate H and 
scale factor a. Since, by definition, one has H = and 
H = a/flo > at the bounce, the null energy condition 
(i.e., , p^ + Pip > 0) is only preserved provided /C > 0, 
hence the choice of positively curved spatial sections. On 
the other hand, as in inflationary scenarios, the strong 
energy condition (namely p^p + Sp^p < 0) is necessarily 
violated. Note that this violation of the strong energy 
condition is nothing but the usual requirement that a > 
in an inflationary stage. 

Following Ref. [2(1, we now switch to a description in 
terms of conformal time as defined above and expand the 
scale factor a, the scalar field (p and its potential V{ip) 
around the bounce, set, for definiteness, to take place at 
r] = 0. To fourth order in rj, we have 



a{r)) 
V{ip) 



3! \m 
1 



1 



(10) 

(11) 
(12) 



r 



where rj^ {agrjo in physical units) defines the duration of 
the bounce and provides a natural time scale. In Eq. ( 10 1, 



the fourth order term with ^ 7^ 0, is written in this way 
in order to emphasize the deviation from the quasi-de 
Sitter solution, namely 0(77) = oqa/I -I- tan^ (v/vo): this 



is discussed in Ref. |2U]. The asymmetry parameter in the 
O (j]'^) term with 6^0 thus corresponds to a deviation 
from a purely symmetric bounce. 

Inserting Eqs. ([lO|, ([TT]) and ([l2]) into Eqs. (|5|, (|6| and 
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Q, and assuming^ (p'g ^ 0, one can relate the scale factor 
parameters to the field values at the bounce as follows, 



3 - T 



(13) 



which gives the value of the scale factor at the bounce, 
and where Vq = V{(po) and T = ^(p'q, as in Ref. PU] . 



1 



1 - T' 

which gives the characteristic duration of the bounce 



(14) 



(3 - T) VT 
3^/2(1 - T)5 ^o' 



(15) 



and 



4T(T- 1) 
5(T-1)2 " 
Vq" 2T(3 - T) 
n 5(T-1)2 



V^y (T-3)2 
Vo) 5(T-1)2 



(16) 



the two parameters describing the deviation from a 
purely de Sitter bounce. Note that in Eqs. ( 15 1 and ( 16 1, 
we have defined a prime, when applied to the potential 
V, to denote a derivative with respect to (p, the function 
being evaluated at the bounce, i.e., Vq = V^{(po)- Note 
that since 77§ > 1, we also have < T < 1. 

As made clear in •20J, if one restricts the analysis to a 
symmetric bounce, one has 6 — such that either ipQ — 0, 
in which case T = 0, or = V^^ = 0. The first case, 
(P'q = 0, corresponds to the exact de Sitter case (ryo = 1): 
the null energy condition is then only marginally violated 
and only the gauge modes of scalar perturbations inter- 
act with the potential term (see Ref. 41J for a detailed 
examination). The second situation, having T ^ 0, de- 
mands that both ifQ and V^^p vanish at the bounce. 

In both cases, one has (p" = ip" [V (ip) ,p' ^H] so that 
requiring a symmetric bounce shrinks the continuous set 
of all possible trajectories to a discrete set, i.e., in a sym- 
metric bounce, the solutions of Friedmann's equations at 
the bounce are a denumerable set. 

Another interesting qualitative result can be obtained 
by expanding the kinetic term ^(^^ around the bounce. 



corresponds to a de 



^ In the case of a symmetric bounce, 
Sitter bounce since 

lim {p^p + Pip) = 

and in the case of an asymmetric bounce, there is singular be- 
haviour in the evolution of the perturbations. 



One gets 



Vk r- I 

T -I- -5-^/2VT(T-3)?7 - 



2T(1 -T) -F2 



V" 

-^T(T-3) 



(17) 



the right-hand side of which must evidently remain pos- 
itive for p real. In the (quasi)symmetric case for which 
one can ignore all Vq' terms, this simplifies to 



1 Fq" T - 1 

2 Vb T-3 



< 0. 



(18) 



Given the bounds on T, this means that Vq" < 0: the po- 
tential must have a convex part, which the field explores 
right at the very moment of the bounce. This restricts the 
possible shapes for the self- interacting potential V (ip) ; 
our specific choice of Eq. (19 1 in {III does indeed fulfill 
such a requirement. Note that in the more generic asym- 
metric case, although we believe that in small deviations 
from the symmetric case, the convex shape of V (p) at 
the time of the bounce remains necessary, no such firm 
restriction can in fact be obtained, and more freedom in 
the choice of V (v?) is expected to be allowed. 

In the case of a general asymmetric bounce, it is worth 
noting that there are two very distinct ways in which ip 
may evolve if the condition ip' ^ imposed in Eqs. (13 1 
to (16) is relaxed. One may obtain such a bounce if (p 
either moves up the potential until (p' = 0, at which point 
it returns towards the value it started from or if it evolves 
from one minimum of the potential to another one. As 
we shall see, the former implies singular behaviour of 
the perturbation equations at the bounce, and so is not 
considered any further in the subsequent analysis while 
in the latter case, one necessarily has ipQ ^ 0, (Pq ^ 
and Vip 0, and the set of solutions is a finite volume 
of phase space dependent on the form of V {(p) only. We 
discuss these points further in the following section. 



III. PHASE SPACE ANALYSIS 

In the previous section, we were able to determine some 
of the properties V {ip) should satisfy for a bounce to oc- 
cur. In this section, further insight into the full dynamics 
is gained by means of a phase space analysis of the evolu- 
tion of both the Hubble parameter and the scalar field in 
the {ip,ip,H) phase volume. Such an analysis is possible 
only once a specific form of the potential V{(p) is given. In 
the remainder of the paper, we consider a spontaneously 
broken symmetry (Mexican-hat) potential of the form 



Vi^) = Vo-^p' + ^p\ 



(19) 
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where /i and A are the mass and self-interaction parame- 
ters respectively, and Vq is the height of V{(f) for ip = 0. 
Such a potential, quite apart from being often used for 
different inflationary models, is theoretically well moti- 
vated for instance as stemming from a grand unified the- 
ories framework, in which case ip is identified with a com- 
ponent of a larger scalar multiplet, the other degrees of 
freedom of which are somehow frozen at the energy scales 
under consideration. 



The potential (19 1 might lead to the formation of do- 
main walls [42, 43J, which, forming in the very early uni- 
verse, would spoil its overall evolution and result in a 
singular crunch [44'. However, if ip belonged to a larger 
multiplet, the symmetry-breaking scheme would hope- 
fully be sufficiently different so as not to produce such 
defects, although even if the decoupling of the various 
degrees of freedom is effective enough that domain walls 
should have formed, it is not absolutely clear whether 
it would have such a dramatic impact. Consider first 
the inflationary case. The maximum temperature, once 
thermalization is taken into account, is found to be quite 
below the Hagedorn temperature, so that the symmetry 
is, in fact, not restored and no defects can be formed at 
all as one proceeds backwards in time ^Sl US]- In the 
bouncing case for which one starts out with a large and 
cold universe, the argument does not stand anymore, as 
the field must be assigned arbitrary values in regions of 
spacetime separated by a distance larger than its corre- 
lation length, in practice its Compton wavelength (note 
that the "causality" argument cannot apply in this 
case, as the horizon can be assumed infinite in the past 
infinity at which one sets the initial conditions). Tak- 
ing into account the effects caused by the evolution of 
the background, it seems probable, though by no means 
guaranteed, that the subsequent contraction would result 
in high decay rates of wall and antiwall configurations, 
possibly leading to a simple and automatic solution of 
the wall puzzle (with tremendous production of parti- 
cles). The remaining wall distribution can, once trans- 
fered into the expanding epoch, be made consistent with 
observations, e.g., the CMB |351 [3D], provided the net- 
work is strongly frustrated |47], or if an initial bias, how- 
ever tiny USl Sni Hn] , is introduced in the potential. 

As it turns out, the phase space analysis is most easily 
performed using the cosmic time variable t. It provides 
some valuable insight, particularly on the space of initial 
conditions that ensure the stability of the background 
spacetime and the occurence of a bounce somewhere in 
its history. As we shall see, a potential of the form (19) 



leads naturally from a contracting phase to a bounce and 
on to a slow-roll type inflationary phase. In what follows, 
we first discuss the set of possible trajectories restricting 
ourselves to the planes (H, ip) and {ip, ip) and then go on 
to describe trajectories in the full volume {H,ip,ip). To 
build the phase portrait in the {ip,H) plane, we start with 
a set of evenly distributed initial values for ip and H, for a 
fixed value of ip, and construct flow vectors with compo- 
nents H and ip defined using (|6| and alj as illustrated on 



Fig. [T] The JC = +1 region is bounded above and below 

by the IC — solution H — ±^ | {ip^ + V), denoted by 

the two parabolic thick curves in Fig. [T] beyond which the 
trajectories describe the /C = — 1 case. A simplified view 
of the full dynamics, as provided by the three slices of 
the figure, readily indicates that a closed universe under- 
going a period of slow-roll type inflation has a nonzero 
probability of having crossed the H = plane in the 
past and therefore undergone a period of slow-roll type 
contraction prior to the bounce, as indicated by the two 
central trajectories (denoted by the dashed curves) in the 
right panel. Projections of typical bouncing trajectories 
in the (H, ip) plane are shown in Fig. |2]for various initial 
conditions. As we shall see in the subsequent analysis, 
there exists, in the full 3D phase space, two saddle points 
located at 77 ~ ±3 and ip = ip :^ 0. In the 2D slices of 
Fig. [l] these saddle points reduce to a repulsor and an 
attractor respectively, not only in the p — slice but also 
in the two other ones, for which the locations of the crit- 
ical points of course shift away from H ~ ±3 to ~ ±1.75 
(center panel) and then collapse into a single point (left 
panel). Given the existence of these two saddle points, 
some trajectories describe a situation in which a universe 
undergoing slow-roll contraction naturally bounces and 
then experiences slow-roll inflation in the usual way. Go- 
ing from right to left in the figure, in the direction of 
increasing ip, one finds that the stable region in which a 
bounce occurs shrinks and the left-right symmetry of the 
right panel is lost. When the field reaches the bottom 
of V {ip), Vnuni the region has shrunk to a point. Some 
other trajectories, on the other hand illustrate a singular 
evolution for which H grows increasingly negative and ip 
is driven to oo. 



The phase portrait in the {ip,ip) plane can be con- 
structed in a similar way. To do so, we fix H and set 
up ip and ip initial values with a flow given by the ip 
values as the vector component in the ip direction, and 
use Eq. ([T]) to define components in the ip direction (see 
Fig. |3] and Fig. |4|. On the top and bottom left panels 
of Fig. [3] -ff < 0, and if p is perturbed away from Vmin, 
both p and ip will be driven to large values. The top and 
bottom right panels show that for H > 0, ip ^ ip^ ~ 2 
(</?zp are the values of ip for which V (ip) is at its mini- 
mum Vmin) and ^ 0, as expected. Note also that as 
— > 0*, trajectories are increasingly circular and that 
in each 2D slice, there exists spiral points at p = ±2 and 
ip = 0. Anticipating on the forthcoming results of the 
3D phase space analysis, one notes that by arranging the 
initial conditions for H and ip appropriately near Knin, 
ip can be pushed towards zero, the region in which tra- 
jectories are nonsingular and undergo a bounce will be 
approached. 

Turning to the full {ip,ip,}i) volume (see Fig. [s] for an 
illustration of what a trajectory might look like) , we may 
now determine both the location of the critical points and 
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FIG. 1: Phase portrait in the {tp,H) plane for /i = 3/2, A = 3/i^/2Vb and fixed values of equal to —2, —1.25, in the 
left, center and right panels respectively. The flow vectors are located at equally spaced values of tp and H, are normalized to 
magnitude unity and their direction is determined using ([6| and |7|. In each panel the isoclines are indicated by thin lines and 

the thick parabolic curves are the = solutions, i.e., H = ±J\ (<^2 ^ yy ^he region in between the curves corresponds to 
K. = 1, while the two remaining regions, one of which is above the top thick curve and the other below the bottom thick curve, 
are K, = —1 regions. The trajectories displayed as thick dashed lines correspond to solutions of the coupled system Q and ([7|, 
some of which connect a slow-roll [ip <^1), H <Q region across H = onto a slow-roll {tp <^ 1), H > region, while others are 
either singular or purely inflationary solutions in the region IC — —1. In the center panel increasingly more trajectories have 
singular behavior as compared to the panel on the right because the region containing nonsingular bouncing solutions shrinks 
as one moves away from (p = 0, i.e., away from the top of the potential V (ip). In the left panel, which corresponds to the 
bottom of the potential V {(p), this region has shrunk down to a single point. 



their stability. The full differential system reads 



dip 
dt 

dV 



A 



dH 



+ 4, - ^' 



(20) 
(21) 
(22) 



The critical (fixed) points are determined by setting 
Eqs. (20 1 to ( p2| ) to and solving for ip, ip, and H. One 
finds 



ip = ± 



ip^± 



(^ = 0, = 0, H = ± 
0, ff=-± 
= 0, H = T 



if : 







3 ' 




(2FoA - 




6A 


{2VoX - 


3/i4) 



6A 



Linearizing Eqs. (20 1 to (22 1 around these six critical 



points and determining the corresponding set of eigen- 
values provides a measure of stability. We first find that 

^0,0, \/Vo/3^ is a saddle point SPl with two negative 



eigenvalues crj^^^-', a-^^^'' and one positive eigenvalue 
a^^^\ given respectively by 



fSPl) _ 



= -2 



3 ' 



(23) 
(24) 



Similarly, ^0, 0, — -\/Vb/3^ is another saddle point SP2, 
this time with two positive eigenvalues crj^^^^ , a^^'^^ and 



one negative eigenvalue a^^'^\ namely 



(SP2) 



(SP2) 



= 2 



3 

Vo 



3 ' 



^(v/3l^±v/3K)+V) 



-2,3 2 

The four remaining critical points have eigenvalues 
(71 = -2Ho, 



CT2,3 



(25) 
(26) 

(27) 
(28) 



where we have set Hq = v/(Vo/3 - /i4/2A). One sees 
that 9H^ - 8/.t2 > follows if A > 9/// [2 (3Vb - 8^^)] , 
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FIG. 2: Trajectories projected in the (<^, H') plane for 
/i = 3/2, A = 3/x'*/2Vb and for four different sets of 
initial conditions on {H,^p,ip), namely (0,0,1/10) [dot- 
ted], (7/10,0,1/10) [short-dashed], (^S/S^Vb/S, 0, 1/5^ [long- 
dashed] and (^-3/5^Vo/3,0, 1/5^ [full]. Note that trajec- 
tories are practically undistinguishable on either side of the 
short bouncing phase during which they cross the H = axis. 



SO that CT2,3 < (resp. > 0) when Hq > (resp. < 0), 
i.e., they are either an attractor or a repulsor. 

On the other hand, 9Hq — < implies that (T2,3 
are complex valued. In this regime, they are either 
asymptotically stable (ASpP) or plainly unstable spiral 
points. The special situation A = 3/i^/2Vb corresponds 
to V {(fi^) = V {(p~) = 0, i.e., a vanishing cosmological 
constant A = 0; the four distinct critical points then re- 
duce to two simply stable spiral points with H = 0. In 
the following, largely for the sake of simplicity and since 
most of the physics we investigate takes place at the top 
of the potential, we set A = 0. 

As already suggested by Figs[l]and|3] bouncing trajec- 
tories necessarily experience both a slow-roll inflationary 
phase, in the neighbourhood of the saddle point SPl and 
a slow-roll contracting phase in the neighbourhood of the 
saddle point SP2. During inflation, the Hubble parame- 
ter reaches its maximum value a/Vq/S (~ 3 in the flgures) 
and then decreases along a ridge while (p and ip will spi- 
ral into one of the ASpP. Evolving the system backward 
in time, one flnds that after crossing the H = plane 
the universe undergoes a slow-roll contracting phase with 
H = — \/Vb/3 (~ —3 in the flgures) followed by a spiral 
towards one of the unstable spiral points. The trajectory 
of Fig.lslis easily understood from combining trajectories 
in Figs^and|3] Indeed, begining with appropriately flne 
tuned initial conditions at Vmin with iy9 ~ — 2, if and 
H < 0, the system can naturally follow a trajectory that 
leads it towards ip = and ip — (going from left to right 



in Fig.[T|, H will therefore become increasingly negative. 
This corresponds to the spiraling motion and the evolu- 
tion along a ridge in Fig. |5] This evolution drives the 
system to SP2, where slow-roll contraction occurs. The 
system then naturally undergoes a bounce, reaches SPl, 
evolves out along a ridge (top right panel of Fig. [s]) , as 
decreases towards smaller positive values. The spiralling 
into the ASpP then ensues. 

At this stage of the analysis, the question of initial 
background conditions naturally arises. Up to now, we 
have identifled a flnite phase space volume in which all 
trajectories bounce, two saddle points and four spiral 
points. Working in this very same volume, one may fur- 
ther deflne the H = surface on which lies a closed 
H — curve. One can then identify a region on the 
portion of the surface bounded by the closed H — 
curve through which the dynamical system can escape 
the volume within and lead to a singular universe. One 
flnds that there is a sizeable region on the H — sur- 
face for which H < 0. Combining volume and surface 
information, one concludes that initial conditions can be 
taken safely in the bulk only, or near the surface but 
in the region for which H > 0. The critical points 

(^0, 0, ±-y/Vb/3^ lie on both the H = surface and the 

H = curve; so do (±y4TVA?,0,0). G iven that the 

latter are spiral points, it turns out to be extremely dif- 
flcult (at least numerically) to choose initial conditions 
leading to trajectories inside the favored volume: ini- 
tial conditions taken in the neighborhood of these critical 
points but only very slightly away have a high chance of 
leading to trajectories that escape the volume and even- 
tually become singular. 

Let us now comment further on the required amount 
of flne-tuning necessary for the bounce to take place. 
Among the set Vtotai of all possible solutions, the sub- 
set that bounce, Vb C Vtotai say, is either of dimension 
-Db — 3 or less (Db < 3) thus leading to a fractal struc- 
ture (see Refs. [H] |52l [53] for such considerations). In 
the latter case, the set of acceptable solutions is of zero 
measure, so one could argue, in the absence of a spe- 
ciflc mechanism that would impose precisely those, that 
a bouncing phase could not have taken place. More work 
is needed to clarify this point. However, in the (more 
plausible?) case = 3, the initial conditions needed to 
initiate a bouncing trajectory can be more easily driven 
to the allowed region (see however Ref. ^^j). Not know- 
ing what happens before the phases we describe here, it is 
essentially impossible to conclude on the amount of the- 
oretical flne-tuning: the situation is akin to demanding 
vanishing spatial curvature /C — 0, i.e., a point (measure 
zero) on the line of possibilities or, invoking inflation, a 
small but flnite region leading to an observationally small 
spatial curvature. 

One point can however be discussed explicitly, and it 
concerns the required level of numerical precision needed 
to reach a bouncing solution starting from a large con- 
tracting universe. Let us assume that we begin the calcu- 
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FIG. 3: Phase portrait in the {ip,ip) plane for fi = 3/2, A — 3/i'*/2Vo and for fixed values of the Hubble parameter H, equal 
to — \/Vb/3 ~ —3, \/Vo/3 ~ 3 in the top left and right panels respectively and —0.1, 0.1 in the bottom left and right panels 
respectively. When H < { left-hand side), the system is naturally pushed away from ip = —2 towards large values of both ip 
and ip. The opposite is true when H > { right-hand side). Combined with Fig. [l] it is clear that initial conditions taken at 
the bottom of the potential, i.e., ifi ~ —2, ifi and H < 0, have to be appropriately (fine) tuned in order for the system to 
stay nonsingular (i.e., ip and (p remaining small) and to evolve towards the region in which a bounce ensues. When H > (top 
and bottom right) the system remains stable and naturally tends towards the bottom of the potential where p = 2. 




-1.0 -0.5 0.0 0.5 1.0 



FIG. 4: Trajectories projected in the {p, p) plane with the same sets of parameters as in FIG.g 
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lation at some time tini say, where we want to impose ini- 
tial conditions. Writing the Friedmann constraint equa- 
tion at tini yields 

with A'ini the number of e-folds we start with. It is clear 
that the denominator is of order e^"" . In other words, the 
larger the universe prior to contraction, the more numer- 
ical accuracy is necessary at iini to satisfy the constraint 
equation. And this is even before one asks the question of 
the Lyapunov stability of the dynamical system. Again, 
all these considerations do strongly depend on the prior 
mechanism of precooling |54j . which might in fact have 
the ability to naturally lead to the nonsingular region. 

IV. SCALAR PERTURBATIONS 




We now consider a typical bouncing solution, such as 
the one depicted on Fig. [6] where the solid and dashed 
lines represent the time evolution of H and a respectively, 
from a slow-roll contracting phase across a bounce to a 
slow-roll inflationary phase. 

Let us now move on to metric perturbations, and con- 
sider the gauge invariant gravitational potential $ in 
the presence of density perturbations Sf. Working in 
conformal time rj for convenience, and in longitudinal 
gauge [55J ^BJ, the scalar part of the perturbed metric 
reads 

ds^ = (r/) [- (1 + 2$) dri^ + {I ^ 2$) j.jdx'dx'] . 

m 

where jij is the background metric of the spatial sections. 
Einstein's equations, to first order, imply 



$"+2 n - 



fc2 - 4/C + 2 W - W 



$ = 0, 
(31) 

where k — ^ n(n + 2) is the comoving wavenumber (the 
eigenvalue of the Laplace-Beltrami operator), while H. 
denotes the conformal time counterpart of i?, namely 
7i = a! I a — aH. 

Using the generalized form of the Mukhanov-Sasaki 
variable Q [SZlEHlEn] 



2/C 



(32) 



and further defining Q — aQ, Eq. (31 1 can be written in 
the form 



z x" z'x' 



g"+( ^--)q'+( ^ + — - — )g = o, (33) 



where 



aip' 



FIG. 5: Four trajectories joining a slow-roll contracting phase 
to a slow-roll expanding phase in the full three-dimensional 
phase space. The sets of parameters used are the same as in 
previous figures. These trajectories are well understood from 
combining the 2D trajectories of Fig. [T] and Fig. [3] For a 
scalar field ip at the bottom of its potential in a contracting 
phase, one has H < 0, if ~ —2 and if ^ 0. Since H < 0, if 
disturbed away from its minimum by a fluctuation, ip will be 
driven farther away from its vacuum expectation value (see 
the bottom left panel of Fig. [sjand also Fig. |4|. For appro- 
priately (fine) tuned initial conditions, the field will be driven 
towards ip — and (p <^ 1 such that (see Fig. [l] going from 
left to right and also Fig. [2| H will grow more negative un- 
ti\ ip ~ 0, ip <^ 1, where a slow-roll contracting phase (with 
p ~ 0, ip <S: 1 a.nd H ~ -y/Vo/3 ~ 3) wiU ensue. The 
trajectory will then almost necessarily follow a bouncing tra- 
jectory and reach a slow-roll inflationary phase followed by a 
spiral down in the other minimum of V (ip) (i.e., it will follow 
an evolution identical to but reversed w.r.t. the contracting 
evolution just described). Note that, although it is certainly 
a possibility, here we do not consider the case in which the 
field turns around because as we will see, the perturbation 
equations are singular for ip — 0. 



and 



z = x 



fc2 + /C 6 



Note, incidentally, that Eq. (33 1 reduces to the usual form 



when the spatial sections are fiat. 



(34) 



It is obvious from Eqs. (^31 1 and (34 1 that the equa- 



tion of motion for <i> is well behaved through the bounce. 
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even though H goes through zero, while that for Q is 
well behaved far from the bounce in the so-called reheat- 
ing (precooling) phase, that follows (precedes) the infla- 
tionary (contracting) phase in which Lp' oscillates around 
zero. Eq. (33 1, however, appears well behaved in neither 



one of these two regimes. 

In fact, Eq. (34) can be made more compact by taking 
zq = aLp' /TL, in which case one obtains the well-known 
form of the equation for the perturbations as that of a 
parametric oscillator, 



0" 



I0_ 



Q-0. 



(35) 



This well-known form is indeed simpler than (34), but 
the potential term is of course still singular when Ti = 
and Eq. [35] is therefore not useful through the bounce. 
The inverse set of remarks can be made of the often used 
variable u, which is related to <i> by 



where 



3n 

2^ 



1/2 



U, 



1 - 



3/C 



1/2 



a \Pf + Pip J V Pipa^ 
leading to an equation of motion for u which reads 



(37) 



fc^- — -3/C(l 



u = 0, 



(38) 



where the potential term is U{t^) = 9" /O + 3/C (l - c^) 
and given explicitely by 



If 



(39) 



Note that this potential grows as a? away from the 
bounce and is not well-defined in the oscillatory preheat- 
ing (cooling) phase that follows (precedes) the bounce. 
The virtue of using the variable u in fact lies in the fact 
that Eq. (38) is well behaved at the bounce, as opposed 
to Eq. ( 35 ) . The form of Eq. ( 38 1 is therefore particu- 



larly well suited to understand the effect of the potential 
term [/(ry) in the neighborhood of the bounce. Indeed, 
substituting Eqs. (14), (15) and (16) into the asymmet- 



ric version of Eq. (48) of Ref. |20). the amplitude of the 
potential at the time of the bounce reads 



f/b = (5-T) + ^(3-T)+(^) (40) 



(T-3)^ 



so that taking the limit T ^ leads to ip ^ Q and (5^0 
and approaches a symmetric de Sitter bounce. Note that 
we have switched notation: what was labelled using the 
subscript 'o' in § [ll| and denoted quantities evaluated at 
the bounce now bares the subscript 'b' in order to dis- 
tinguish these quantities from Vq which now appears ex- 
plicitly in the following expressions as the constant term 



of the potential V {^p) 
has 



V{(p = Q). Using Eq. (IT9| one 



3 (mV^ - 4V-o) 
12^2^ - 2\ip^ 



and 



b 



12 



5^2 
6A 



6Vo - 5Ai2 



(41) 



(42) 



such that if the scalar field goes to either one of the limits 



±1 



6/i2 ^ 2 (V - e^oA) 



1/2 



± 



A 



(43) 



the potential exhibits a large central peak. However, for 
the theory to make any sense at all, the potential should 
be positive-definite. This is ensured provided one de- 



mands Knin > 0, where V,n 



Vb - 3/i7(2A) is the 



(36) 

minimum value of V , attained for tp — > yJ&Jj^JX. This 



means that A > i^J.'^/{2Vo). In Eq. (W3| however, the 



argument of the square root requires the opposite condi- 
tion A < 3/x^/(2Vo) to hold in order to keep tp real. The 
abovementioned limit can therefore not be approached 
without violating one or the other condition: the quanti- 
ties V{^/V\j and IV\^ must therefore remain small. The 
potential Uh, having no large central peak, therefore does 
not affect the perturbations much as they pass through 
the bounce itself. The question of the wings (located at 
^ ±0.75 in Fig. [7), and in particular whether there is any 
way of showing they remain small, is left unanswered. 

It should be noted that the variable Q is equivalent 
to the variable often denoted by v in the literature (see, 
e.g., Ref. |56jV It arises naturally when expanding the 
Einstein action up to second order in the perturbations, 
giving the action of a scalar field in fiat space having 
a varying mass (i.e., a pump term) caused by gravity. 
It allows an easy definition of the initial quantum state. 
Assuming adiabatic vacuum in the low-curvature regime 
for which k'^ » z" / z in Eq. (35), one sets 



^init oc 



ikrj 



2k 



(44) 



We have already seen in Section |III| that the initial con- 
dition for the background is quite difficult to achieve 
in a natural way. We now encounter another weakness 
of curvature-based bouncing models such as this one, 
namely that there does not seem to be a generic way of 
setting up appropriate initial conditions for the perturba- 
tions. This is due to the fact, revealed by Eq. (34), that. 



because the potential U never drops to negligibly small 
values far from the bounce, the condition fc^ ^ z" / z is 
never satisfied because a^y^^ dominates. 

What, then, should we choose as initial conditions for 
the perturbations? 

In slow-roll infiation, <i> = const, in the regime k <C aH , 
such that one has <i> = 0. Fixing an initial condition in 
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FIG. 6: Number of e-folds, log a, (dashed-dotted) and Hub- 
ble parameter H (full) as functions of time, for a symmetric 
bounce, with the Mexican-hat potential V [ip) of Eq. (19 1, 
with parameters Vo = 1, = 1 and A = 3/i*/2Vo. In order 
to obtain this solution, initial conditions were taken at the 
bounce, with (/jq = 1/10. Note that these values were chosen 
such that the Hubble parameter remains almost constant in 
both the contracting and expanding phases for ~ 70 e-folds. 



this regime seems hardly motivated indeed, as it would 
be equivalent to assuming an initial spectrum, contrary 
to the basic original idea of computing this spectrum out 
of a natural prescription. However, the model we are con- 
sidering here cannot be made complete, in the sense that 
a contracting epoch should be connected to the bounce 
through a yet undescribed precooling phase [S^ . Let us 
then press on and assume the usual form 



$ oc C{k) 



1 - 



H 



, it) dt 



(45) 



in which the lower bound ti, on the integration is 
unknown and we have parameterized our ignorance 
through the unknown overall multiplicative function of 
the wavenumber k, C(k). 

As first recognized in [3Dj, Eq. (31 ) admits a first inte- 
gral given by 



Cbst 



' ■ •? 



1 



such that [5T] 



Cbst oc 



It5S 



(46) 
(47) 



where the nonadiabatic pressure perturbation is given in 
terms of the entropy perturbation 5S through 



1 



o'tSS = 6p^ - cl6p^ = (1 - cl) (3/C - e) (48) 



FIG. 7: The potential U (solid line), and its approximation 
(dashed-dotted) in the slow-roll region, together with i^bst 
(dotted) as functions of time, obtained from the evolution 
of u using (381 in the background of Fig. |6] with n = 20. 
We should emphasize that, even though this is not immedi- 
ately obvious from the graph, ^bst is not conserved across 
the bounce. 



with defined by the adiabatic variation of ptp and 

as 



(49) 



and is the square of the sound velocity, i.e., the velocity 
at which perturbations propagate. 

For a single scalar field in flat space, 5S oc k^^, as a re- 
sult of which the entropy perturbation is essentially neg- 
ligible with respect to the adiabatic ones, and C, vanishes 
as long as fc <C a. The curvature perturbation C thus re- 
mains constant on super-Hubble scales. In the inflation- 
ary scenario, one uses the variable v and the dynamical 
equation for the perturbations can be solved directly in 
terms of Bessel functions of the first kind of order a where 
a is a constant that can be expressed in terms of the 
slow-roll parameters e and (5, namely a ~ —3/2 — 2e + 5. 
Using the asymptotic behaviour of the Bessel function 
and performing two consecutive matchings with the vac- 
uum solution of Eq. (44 1 for krj — s- —oo on the one hand, 
and with a solution of the type (45 1 for fcr] ^ on the 
other, one obtains 



= fc3|C|2 ^ k^A^{k)k 



2a+l 



(50) 



where A{k) is the power law k dependence of the ini- 
tial state. One sees that the standard slow-roll spectrum 
^-4e+25 g|^g]-Qg directly from the Minkowski vacuum, in 
which A{k) ~ /e-i/2. 

In the bounce scenario, the situation is slightly dif- 
ferent, as, in particular, no specific prescription for the 
initial state exists. Given an appropriate potential U, 
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one may however similarly perform several consecutive 
matchings. We have already shown that the potential 
for u at the bounce is negligible. This means that the 
solution in the neighborhood of the bounce can be taken 
as a linear combination of oscillatory functions, with un- 
known /c-dependent coefficients that depend on the ini- 
tial state. In the slow-roll regions on either sides of the 
bounce, where curvature can be neglected, the solution 
is a linear combination of Bessel functions while very far 



from the bounce, where k <^ aH , it is given by Eq. (45 1 



In practice, it turns out to be sufficient to perform only 
two matchings. The first joins the initial conditions (for 
which a reasonable prescription has yet to be given) to 
a solution in terms of Bessel functions, while the sec- 
ond joins Bessel functions solutions having arguments 
k\rj — ry_| (respectively k\'q — rj+\) and order (resp. 
(3+) where rj± and a± are defined using the approximate 
form of the potential far from the bounce. In the slow- 
roll regions, for which the curvature term is unimportant, 
being suppressed by a factor a?, the potential for u can 
be expressed as 

U a^H"^ {2e-5 + 2e^ + 5"^ - ie5 - C) (51) 

where e, 5 and ^ have their usual definitions in terms of 
the derivatives of the potential w.r.t. the field [HH 
In such a case, the potential term before {U-) and after 
(C/+) the bounce reads 



a± 



(52) 



where r]_ and ?7+ are the asymptotic limits to the left 
and right of the bounce respectively. In de Sitter ri± = 
±7r/2. The indices a± are obtained from the slow-roll 
parameters e±, S± and ^± as 



4a| = 1 + 4 (2e± - d± + 2e| + - 3e±5± - £_±) 



(53) 



which, for small values of e±, 6± and ^± gives 



(2e± -S± + 2el + 4 - 3e±'5± 



(54) 



For a symmetric potential V{ip), provided sustained slow- 
rolling contracting and expanding phases are achieved, 
a_ and a+ are equal, because Fig.|5]is symmetric w.r.t. 
the {if, if) plane. This is clearly true in particular in the 
case of a symmetric bounce, but most likely also even in 
most asymmetric bounces. 

In order to set up sensible initial conditions away from 
the bounce, i.e., where k'^ <C U, one has, using Eqs. (36 1 



and (45) 



C{k) 



1 - 



n 



(55) 



Here, C{k) cx k'', and both 7 and 77* are constants that 
depend on a preceding precooling phase. In the de Sitter 
limit, one obtains the approximate initial condition 



Mi cx C(fc) 



Note also that far from the bounce, 



,BST 



n 

aip' 



(56) 



2^^, (57) 



so that Q oc C{k), a mere function of the wavenumber 
fc, i.e., a constant in time, the k dependence of which is 
the same as that for u. The matching is performed from 
a slow-rolling contracting phase through the bounce into 
an expanding slow-rolling phase, i.e., from 



-(??) - ^k{T^-T^^){AiJ^_ [(/c(?7-?7_)]+A2J-o_ [fc(ry -??-)]}, 



where 



and 



to 



where 



2-^+"-7r[fc (?7.-r/-)]- 



-l/2-Q_ 



[ui (2a_ - 1) + 2u[ {rji, - ?7_)] esc (7ra_) 



-1/2-1-0 



r(i-a_) 



[u, (2a_ + 1)- 2u[ (7,, - r?_)] T (-a_) , 



, (77) = ^k {v-V+) [Bi (fc) Jc.^ (fc (7? - 77+) ) + i32 (k) J_„^ (fc (ry - 77+)] 



Bi (k) 



Aicos (77. - ?7+) + - (a_ - a+) 



A2 cos 



k (r/_ - 77+) - - {a- + a+)J j 



(58) 

(59) 

(60) 
(61) 
(62) 
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and 



B2 (k) = CSC (yTTa^) cos 
The mixing matrix is thus given by 



k iv- -11+)+ 2 + "+) 



A2 cos 



TT ] 1 



(63) 



CSC {^^a-^.) 



cos [fc (?7_ - 77+) + I (a_ - «+)] cos [fc (?7_ - 77+) - | («_ + «+)] \ / Ai 

cos [fc (77_ - 77+) + I (a_ + «+)] cos [fc (77_ - 77+) - | (a_ - 0:+)] J \ A2 

I 



This result differs from the one obtained in [20, in two 
ways. First, the characteristic time over which the per- 
turbations may be altered by the bounce: in [20], it 
was extremely small, and the effect was mainly concen- 
trated at the bounce itself ^, whereas in the present case, 
taking into account the quasi-de Sitter contraction and 
expansion phases that take place before and after the 
bounce, the actual characteristic time scale turns out to 
be rj- — rj^. The second crucial difference is that the 
bounce itself leaves the perturbations unchanged, since 
the potential is negligible during this stage of cosmologi- 
cal evolution. Our result however agrees with Ref. [20 in 
the sense of Eq. (70) of the abovementionned reference, 
i.e., in the immediate neighborhood the bounce, one gets 
the identity, together with some possible overall rescal- 
ing of the total amplitude. We however expect the result 
presented here to hold more generically but we remind 
the reader that as already mentioned, the question of the 
amplitude of the wings of the potential is left unanswered. 

Let us turn to the power spectrum. For fc(|77— 77_|_|) ^ 1, 
the first term is negligible, and the result can be approx- 
imated by the term in 82- In addition, A2 also turns out 
to be negligible. One thus has 



Vc (X fcl+27-25o 



-2So 



COS" yt (77_ - 77+) -I- - («_ + Q+) 

(64) 

where k'^ is the initial fc dependence of Ui and u[, i.e., 
that of C{k). The power-law coefficients 5a± are given 

by 



(65) 



From the above, we find that as 1^ gets vanishingly small, 
then TT — (7;- — 7;_|_) — > TT^, and both ba± — > as /i is 
decreased. In the special case of 7 = —1/2 (in which case 
one starts out with an almost scale-invariant spectrum 
for Vc^ far in the past prior to the bounce), the resulting 
power spectrum takes the form 



oc 



-25a_-2i5Qa 



COS 



k{ri- - V+) + 2 ("-+"+) 

TT 

(tt + Sr]) H — {6a- + 



(66) 




FIG. 8: Expected spectrum of the curvature perturbation 
Cbst. The theoretical prediction of Eq. ( |66[ ) is the sohd Une, 
while its approximation of Eq. (68 1 is the dashed line. The 



crosses are t he c orresponding numerical results obtained by 
evolving Eq. ( 38 1 for wavenumbers n in the range n £ [60 — 
2000]. The power law term alone is shown in the dashed- 
dotted line. Here, the parameters for the potential V{i.p) of 
1/5, if 
and Srj ■ 



Eq. ( 12 1 are ^ ; 
5a+ — Sa^ ~ 10 



1/30; these parameters imply 
10"" 



where we have set, for notational simplicity, 
Sr] = n- (77_ -77+). 



(67) 



Taking fc = y/n{n + 2) 
approximate form 



1 — l/2n, one obtains the 



which we compare, on Fig. |8] with Eq. ([66| and with the 
actual full numerical solution. It is clearly a satisfying 
approximation as the points and curves are almost undis- 
tinguishable. Not only does the bounce modify the power 
law behavior of but it also adds to it a multiplicative 
oscillatory term at leading order. In first approximation, 
one further has 



Sr] 



1 



^a^Vo ~ 2/C 



- 1 



(69) 



^ A (5 function approximation for the potential was even suggested. 



This highlights the fact that this bouncing scenario is 
fully characterized by the usual slow-roll parameters en- 
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coded in 5a± and by the scale factor at the bounce, ag, 
that appears in the expression for ^77. The above oscillat- 
ing spectrum can be rephrased, once the inflation epoch 
has taken place and ended, into a more generic form 



(70) 



where now the wavenumber fcph is written in units of in- 
verse length, the nominal scale fc^ can be fixed at the 
given scale of 100 Mpc and (p is some phase. One is thus 
left with one arbitrary parameter, namely the oscillation 
frequency uj. For very small frequencies, i.e., oj <C 1, the 
cosine term remains constant and this is nothing but a 
simple power-law spectrum, and no new information can 
be derived. On the other hand, in the large frequency 
limit, w ^ 1, the resulting spectrum, once evolved in the 
radiation and matter dominated eras and expressed in 
terms of CMB multipoles, is not immediately ruled out, 
contrary to what might have been expected first hand. 
Fig. |9] shows that for the present best set of data from 
first year |Mj or third year WMAP [BS] , a very reasonable 
fit can be obtained [66j El- Clearly, with a primordial 
spectrum as the one exhibited on Fig. |8] one might have 
naively thought it impossible to recover a temperature 
perturbation spectrum such as shown on Fig. |9] It turns 
out, however, that provided the frequency uj is sufficiently 
large the spherical Bessel smoothing functions used to de- 
rive the latter from the former j68j allow for a fit which 
is neither compellingly wrong, nor any better that the 
usual power-law spectrum. Similarly, in the observable 
range of scales where only linear interactions play a role 
in the eras following inflation, and for which the transfer 
function from the primordial power spectrum is a sim- 
ple function, oscillations in the density power spectrum 
can nonetheless be expected to be smoothed out by a 
convolution with a windowing function linked to survey 
size in observations. Constraining the parameters of this 
model using large-scale galaxy surveys and CMB data 
will be discussed in future work. In particular, relating 
the parameters and A in Eq. (70 1 to those of our 



spectrum ( 68 ) , and to those appearing in the potential 



(19 1, namely Vq, fj,, and A will be needed to go any fur- 



ther and place meaningful constraints on this particular 
type of model. 



V. CONCLUSIONS 

Using the simplest possible theoretical framework 
in which a bouncing Friedmann-Lemaitre-Robertson- 
Walker scale factor is possible, namely GR, positively 
curved spatial sections and a scalar field in a self- 
interaction potential, we investigated the set of initial 
conditions under which a bounce might occur. For the 
bounce to take place, we assumed a potential of the form 
( 19 1, with a constant term Vq, a symmetry-breaking mass 



parameter fi and a self-interaction coupling constant A 
and without loss of generality in this work, A = 0. This 
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FIG. 9: Using the spectrum of Fig. [S] is not obviously in- 
compatible with the data: WMAP data (1 and 3 years, re- 
spectively, squares and circles) are superimposed on the pre- 
dicted (full line) temperature fluctuations stemming from 1 70 1 
for fixed fc*, lj = 2000 and 1000 for the upper and lower 
plots respectively and ip was set to zero for simplicity. Both 
the amplitude and frequency of oscillations are affected by 
varying the parameters of the model. Figure courtesy of C. 
Ringeval [69j . 



potential belongs to the more general, low-energy, phe- 
nomenological class of potentials in which all powers of 
the field could be present, including, e.g., cubic terms. 
Such asymmetric terms are however not privileged, and 
we have found that in any case the main results of this pa- 
per are roughly independent of those extra terms. This, 
we interpret, is mostly due to the fact that all the physics 
of the bounce as well as that of the slow-roll contraction 
and expansion phases occurs near the top of the poten- 
tial, the shape of which is hardly dependent on anything 
but the parameters we have used. The background we 
looked at here is therefore, as far as we can tell, generic 
enough that one be in a position to draw general conclu- 
sions on the evolution of the perturbations through the 
bounce, which was the main concern of the second part 
of this work. 

Our first definite conclusion concerns the background 
and is provided by the phase space analysis. It is found 
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that, given positively curved spatial sections, the poten- 
tial V {(fi) as given by Eq. (19), and assuming the oc- 



curence of an inflationary phase some time in the history 
of the universe, a bounce can indeed have occured in the 
past with nonzero probability, preceded by a slow-roll 
contracting phase. Note that the form of the potential 
we have chosen is the simplest one in which a bounce 
can happen and with which ip stays bounded from below. 
The flipside of this result is that given "natural" initial 
conditions in the far past, i.e., a large but contracting 
universe with a scalar field slightly disturbed away from 
the minimum of its potential, trajectories for the evolu- 
tion arc highly unstable and fine-tuning is a necessity if 
one is to head towards a bounce and avoid singular be- 
haviour. The shape of the potential must also be ajusted, 
not unlike what is required in inflation, to ensure an in- 
flationary phase lasting at least 60-70 e-folds. At this 
stage, one should also mention that a slowly contracting 
precooling phase, during which, say, a radiation domi- 
nated universe moves towards a scalar field dominated 
universe, by some kind of condensation process, is a nec- 
essary step that should occur prior to the bouncing mech- 
anism that we have described in this work. The question 
of whether, after such a precooling phase, the scalar field 
dominated /C = 1 universe will follow a bouncing rather 
than a singular trajectory is left open, not only in view 
of the first remark made in this paragraph, but also be- 
cause one can choose to constrain the likelihood of such 
nonsingular trajectories by looking for their signatures 
in cosmological data (CMB and large structure forma- 
tion for instance). This naturally leads to the study of 
the perturbations and attempts to draw generic conclu- 
sions on their behaviour across the bounce for a given 
class of models. 

Curvature perturbations, encoded in C (or Cbst the 
/C = 1 case) are often advocated to be conserved, i.e., 
constant, in several of the known cosmological phases, 
such that they can be used as reliable estimators. Such 
a proposal frequently meets with strong disagreement. 
In Ref. To] for instance, it was shown that out of the 
three reasons why C-bst usually conserved, two at least 
are not fulfilled, and the third is presumably meaning- 
less in the cosmological bounce context. In the case at 
hand, and working in the observable wavelength (large) 
range, we found that Cbst indeed not conserved, with 
both its spectrum and amplitude modified. We used a 



matching procedure connecting a Bessel solution for the 
evolution of the variable u prior to the bounce to an- 
other after the occurence of the bounce. Given an ini- 
tial quasiscale-invariant spectrum, the transfer function 
for the perturbations preserves the initial power-law to 
within order of the slow-roll parameters, but adds a lead- 
ing order oscillatory term. This result was confirmed by a 
full numerical evolution of the dynamical equation for u. 
A priori, the form ( [70| might lead to the hasty conclusion 
that this class of models can safely be excluded by CMB 
observations. We found however that in a large part of 
the parameter space in which a bounce is realized (see 
Section III I, the spectrum of CMB multipoles (Fig.|9| is 
merely modified in a way reminiscent of superimposed os- 
cillations, a possibility that has attracted some attention 
in recent years [33^ , and was found to be not inconsis- 
tent with the data. Any oscillations in the matter power 
spectrum or CMB multipoles are expected to be mostly 
smoothed out either by transfer functions or observa- 
tional windowing functions. One is then led to predict, 
as in the completely different trans-Planckian case [66] . 
some amount of superimposed oscillations which is not 
any worse (and might in fact be statistically better) than 
the power-law only inflationary prediction. Constraining 
the parameters of the whole class of effective 4D bouncing 
models in which the potential term for the perturbations 
at the bounce has a negligible effect, using large-scale 
galaxy surveys and CMB data, remains to be done |71j . 
Once this is achieved, a complete model will need to 
be constructed, starting with an initial large contract- 
ing phase, with appropriate (physically justified) initial 
conditions, filled with matter and radiation, and evolving 
in such a way as to connect to the present model via a 
precooling phase in the recent fC— bounce proposal [27]. 
Such a mechanism is presently under investigation [54^ . 
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